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Synopsis 

A comprehensive, semi- quantitative model for the thermodynamics of hydrophobic sol- 
vation is presented. The model is based on a very simple premise suggested by the scaled 
particle theory and treats both solute and solvent molecules as hard spheres. A connec- 
tion between the peculiarly large heat-capacity change for hydrophobic solvation and the 
large temperature dependence of the thermal expansivity of water is found. Analysis reveals 
a possible physical origin for the converging behavior of solvation entropies for a series of 
homologous hydrophobic compounds. The model suggests that the low solubility and the 
large heat-capacity change of hydrophobic solvation stem from two distinct aspects of water 
molecules: the static geometry of the molecule and the dynamic hydrogen bonding network, 
respectively. 
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Introduction 



The scaled particle theory (SPT), a statistical mechanical model for liquids and solutions 
made of hard spherical molecules, is a powerful conceptual framework and a useful com- 
putational tool for molecular solvation. 1-4 Even though the SPT was originally developed 
only for hard spherical molecules, 5 recent studies have shown that it also provides insight 
into the solvation of hydrophobic solutes in associative liquid like water. There are several 
reasons for this wide applicability of SPT: 

1) The SPT is not a statistical mechanical theory based on the first principle. Rather, 
an independently obtained solvent density has to be provided in addition to temperature, 
pressure, and molecular parameters. This is in contrast to a genuine statistical mechanical 
theory of liquid, in which density is derived from temperature, pressure, and molecular 
parameters. Some essential properties of an associative liquid are contained in the density 
data. 5 ' 6 

2) Recent studies have shown that the SPT provides accurate numerical calculations for sol- 
vation free energy of hydrocarbons in water, but does so with less success for calculating the 
entropy and enthalpy. 2 ' 6 This result is consistent with the concept of entropy-enthalpy com- 
pensation, which occurs with the solvent reorganization when introducing a solute molecule 
into a solvent. 6-9 

3) To apply the SPT to an associative liquid like water, one practical difficulty is how to 
determine the hard sphere radii for water and the solute molecules. 2 While these difficulties 
affect quantitative numerical calculation, it does not alter the qualitative physical insight one 
obtains from the SPT. In this paper, our primary interest is in the qualitative relationship 
between different physical quantities, the numerical values are not essential for our analysis. 

Based on these arguments, we try to use SPT as a semi- quantitative model for under- 
standing various experimental results which are essential to the solvation of small organic 
molecules in aqueous solutions. We show that much of the characteristics of hydrophobic 
solvation can be understood in terms of this over-simplified model. Our approach is very 
much in the same spirit of the earlier work of Grunwald. 8 Treating the hydrophobic solutes 
as inert hard spheres, of course, completely neglects the solute-solvent interaction. 3 ' 6 This 
should be kept in mind when applying the result of the present analysis to the solvation 
of hydrocarbons. On the other hand, the dissolution of inert gases like xenon exhibits all 
the characteristics of hydrophobic solvation. Hence, neglecting solute-solvent interaction is 
justified in our present study. 

Some Basic Facts 
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The essential features in the data for dissolution of hydrophobic compounds in aqueous 
solutions are summarized as follows: 

1) Dissolution entropy changes, A5'(T)'s, associated with transferring a series of hydrophobic 
solutes from solid — > water, liquid — > water, and gas — > water form three distinct, respective 
groups. Within each group, AS'(T)'s for different solute species converge to a common AS** 
at a temperature T*. While AS^'s are different for the three groups, T*'s for all three groups 
are approximately equal to 110°C. 10 The two groups of data on the transfer from liquid — > 
water and gas — > water are in agreement with the Trouton's rule. 

In order to compare free energies calculated from SPT with experimentally obtained 
solubility, one has to calculate the solvation free energy from a measurement based on con- 
centration scale in molarity rather than in mole fraction. 5 ' 11 ' 12 The numerical difference 
between the two calculations is —RTln{vi/v2) where V\ and v<i are the molar volumes of the 
solvent systems before and after the transfer, respectively. For the gas phase, v = ksT/p 
where p is the gas pressure (= 1 atm); for water, v =18.08 cm 3 /mol. Q This changing of 
concentration scale increases the free energy values given by Murphy et al. 10 by about 60 J 
per degree per mole 13 (e.g., shift the lines in their Figure 1 upward). 

2) The heat-capacity changes, AC p 's, for all the dissolution reactions are approximately 
temperature independent. 14 

3) There also appears to have a convergence temperature for enthalpy changes in disso- 
lution of organic compounds from gas — > water (see Appendix for more discussion). The 
temperature is not known, but the corresponding AH* is about a few KJ/mol. 

There have been many different ways of presenting experimental data on the thermody- 
namics of dissolution. Here we have summarized several key results from the literature which 
are useful in organizing and relating various data representations. It is important to note 
that many presentations are in fact concerning same experimental measurements. They are 
related by the basic thermodynamic formulae, as shown in the following theorems. 

Baldwin's Theorem 

Based on the observations of Sturtevant 15 and Privalov 14 , R.L. Baldwin gave the following 
theorem in 1986. 16 For a series of homologous reactions, designated by i, with temperature- 
independent AC Pi , the following two statements are equivalent: 

1) At one temperature T\ there are unique constants a and b such that: 

AS i = a + bAC Pi (1) 

2) there is a unique temperature T* at which AS(T*) is unique. That is, T* and ASi(T*) 
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are independent of i. Hence: 



As t (T) = as(t:) + ac p >(t/t;) (2) 

Compare Eq. (1) with (2), we have a = AS(T*) and b = ln(Tl/Tf). From now on, we will 
follow Lee 13 and refer to the ASi versus AC Pi relation in Eq. (1) as a SMPG plot named after 
Sturtevant, Murphy, Privalov, and Gill. 10 ' 15 Baldwin's theorem states that a linear SMPG 
plot is equivalent to having a convergence temperature for ASi, an d the intersection a in 
the SMPG plot is the unique AS'*. Conversely, the existence of a convergence temperature 
indicates a linear SMPG plot at every temperature. 



Lee's Theorem 

Lee's theorem 13 is an application of the mathematical property of a bilinear function: 
if a function f(x,y) is linear as a function of either x or y, as well as a function of their 
product, xy, then there exists a y* at which f(x,y*) is independent of x. To express this by 
equation: 

r i \ r t u c \i b.ad — bc 
f{x,y) = a + bx + cy + dxy = d{x + -){y + -) H — 

so when y = y* = —b/d, f(x,y*) is independent of x. 

The bilinear function has to be a linear function of either variables when the other one is 
at a fixed value. To apply this result to the thermodynamics of dissolution, we identify the 
logarithmic temperature and a certain molecular parameter as the variables x and y, and 
the entropy (or enthalpy or Gibbs free energy) of dissolution as the function f(x,y). Hence, 
at a given temperature, the entropy is a linear function of the molecular parameter X. Lee 
proposed the parameter to be a size measure of the solute molecule. In fact, Lee's proposal 
of 

AS,(T) =a s + b s X l (3) 

is a structural interpretation of the linear SMPG plot given in (1). The subscript "s" here 
stands for entropy. (There is a similar relation for enthahlpy, with its respective and bh- 13 ) 
It is well known that AC p 's are proportional to molecular size of a hydrophobic solute. It 
is also worth noting that a linear relationship due to substituents of a series of functional 
groups is widely observed in organic reactions. 17 

Therefore, the bilinear argument immediately leads to the existence of a convergence 
temperature from the linear entropy relation in Eq. (3). Combining Lee's and Baldwin's 
theorems, it can be shown that the linear entropy relation is sufficient but not necessary for 
generating a convergence temperature. 



The BMDW Theorem 
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This theorem, which was given by Baldwin and Muller, 18 and also independently by 
Doig and Williams, 19 establishes an intrinsic relationship between the three convergence 
temperatures for entropy, enthalpy, and Gibbs free energy. Consider a series of homologous 
reactions designated by i and each has a temperature-independent AC Pi of its own. If two 
out of the three thermodynamics quantities (entropy, enthalpy, and Gibbs free energy) have 
convergence temperatures, then there is a convergence temperature for the third quantity. 
If we denote by T*, T£, and T* these convergence temperatures, we have: 

t; - t* - t;iu{t;/t* s ) = o (4) 

and also a relation for the corresponding thermodynamic quantities: 

AG(T*) = AH(T*) - T;AS(T S *) (5) 

This theorem was used by Baldwin and Muller 18 to explain the intriguing fact of approx- 
imately equal convergence temperatures T* and T£ for a set of protein folding reactions. It is 
discovered that this set of proteins have approximately equal melting temperatures, i.e., T* 
= 331 ± 9 K and AG(T*) 0. Using the same theorem, Doig and Williams 19 have reached a 
similar conclusion. They pointed out that most proteins of the same set have approximately 
equal AG per residue at room temperature. These two arguments are consistent with each 
other if we note that all the proteins in the set have approximately same AC P per residue. 

A Physical Model Based on SPT 

The basic results of SPT are summarized here. It is not necessary for the readers to 
know the technical details of SPT; 5 rather, our model starts with Eqs. (6) and (7) below. 
SPT is a rigorous statistical mechanical theory for liquids made of hard spheres. The theory 
also provides an approximated formula for calculating the free energy of dissolving a hard 
spherical solute from gas phase into a hard spherical solvent, which is equivalent to intro- 
ducing a spherical cavity in the hard sphere liquid. The free energy, which is a function of 
the radius of solute as well as the density and the radius of pure solvent, has three dominant 
terms. The first term is independent of solute size; it is associated with the reduction of 
conformational space of the liquid upon introducing a solute molecule of zero physical size 
into the solvent. With the physical point being present in the middle of the solution, no sol- 
vent molecule can occupy the same point in space. Hence such conformations are no longer 
accessible. This term is a function of the density and the radius of pure solvent only. The 
second term is proportional to the square of the radius of the solute, i.e., its surface area. 
And the third term is proportional to the cubic power of the radius, i.e., its volume, and 
external pressure. It has been repeatedly demonstrated, both theoretically and empirically, 
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that the volume term is negligible for any molecular size cavity, and the dominant effect is 
from the surface term which defines a proportional coefficient called surface tension. For 
hard spherical solvent and solute, SPT gives an expression for the surface tension, and also 
a minor correction term on the free energy due to different curvature of the solute molecule. 
The free energy is given as: 5 



AG , „ m r 
— = -ln + - 

KbT a 1 



3£ 
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1-C 2 \ 1-f 



(6) 



where k B is the Boltzmann constant, T is temperature in Kelvin, a is the radius of solvent 
molecule, £ is the packing density for pure solvent (the volume fraction occupied by the hard 
spheres), and r is the radius of solute molecule. 
Our model starts with a very simple premise: 



AG(T)=5(T,0 + °(T,Or 2 



(7) 



where 8 and a are solvation energies for a point and surface tension, respectively. Thermo- 
dynamics for solvation will depend on the functional form of 5, a, and implicit temperature 
dependence of £. Eq. (7) offers the possibility for connecting thermodynamics of solvation 
with the thermal expansion coefficient of pure solvent (see below). 
From Eq. (7), it is straightforward to obtain: 



AS(T) = -St - - ((Jt + ^£r)r 2 
AH(T) = 5-T5 T - TS^t + (<t - Ta T - Ta^ T )r 2 



(8) 
(9) 



where subscripts "T" and represent partial derivatives with respect to these variables. 

Let's now estimate the magnitudes of the various terms using values for water 3 (SPT 
parameters for water at room temperature are £ = 0.363 and a — r w — 1.38 A.) £t is 
proportional to the thermal expansion coefficient of the pure solvent: £r = — cn£, where a 
for water is very small. At 1 atm, the values of a for water range from —0.064 x 10~ 3 at 
0°C to 0.7 x 10~ 3 at 100°C, and a = 0.257 x 10~ 3 at room temperature (25°C). The thermal 
expansion coefficients for non- associative liquids are around 1 x 10~ 3 at 1 atm and between 
0°C to 100°C, so these values are not very different from that of water. 

Quantitatively, compare Eqs. (6) and (7), we have: 



k B Ta^ 



k. 



'300 x 0.257 x 10~ 3 x 0.363 s 
1 - 0.363 



0.04fc B 



while 



S T = -k B ln(l - = 0.5A; B » \6g T \ 
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where we have used T = 300K. Similarly by simple differentiation, 

^ = kBT ll{l-ff = 20kBT ^ 
Hence \ct^C,t\ ~ 0.56/cs/r 2 , while ot = ?>.YlkBlr 2 w » \(t^t\- 

The Heat Capacity AC P 

The heat capacity can be obtained from either Eq. (8) or Eq. (9): 

AC P = -T(5 TT + 25 T ^ T + 6#(%. + 5^tt) ~ T(a TT + 2a T g T + a^ T + a^ TT )r 2 (10) 

While the values for a are not very different between water and non-associative liquids, 
there is a dramatic difference between a's dependence on temperature. This difference 
contributes a large term to the AC P of solvation in water, in contrast to non-associative 
solvent. In other words, even though a for water is very small, its temperature dependence 
is quite large, in contrary to most organic solvents. Hence we will neglect contribution from 
£t but shall keep the terms with ^ TT . Thus we have: 

AC p = -T^ TT (5 ( + a e 2 ) (11) 

where, according to SPT (i.e., Eqs (6) and (7)), both 5% and are explicitly proportional 
to T. Therefore if £tt is proportional to T~ 2 , then AC P will be approximately temperature 
independent. 

It is clear from the above argument that, for a non- associative solvent, the solvation of 
a inert solute should have very small AC P since non-associative liquid has almost zero £tt- 
For water, 5% = 1.57k B T and o^r 2 = 20.0k B T(r /r w ) 2 . Therefore, Eq. (11) gives a linear 
relationship between AC P and molecular surface area r 2 with almost zero intersection when 
r > r w . Note, however, that for hydrocarbon solutes like propane and isobutane, there are 
significant solute-solvent interactions which contribute to the overall AC P . 6 Therefore, our 
present result is only semi-quantitative and has to be augmented with such interactions when 
applied to real experimental data. 

AS* and AH* 

We now return to Eq. (8). According to Lee's theorem: 13 

AS* = -S T - 6g T (12) 

Note that in order to compare calculation from SPT with experimental measurements, we 
have to obtain the entropy from experimental measurements according to molarity con- 
centration scale. 5 ' 11 When this was done, 13 Lee found that the calculation given by SPT 
compares favorably with experimental results. 
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Similarly, we have: 

AH* = 5 -T5 t -T8£t= -T5£ t (13) 

the second equality is because 5 is a linear function of T (see Eqs. (6) and (7)). Numerically 
—TS^t is about 0.034 ksT, that is, 2.5 kJ/mol. It should be noted that Eq. (13) neglects 
contribution from solute-solvent interaction. It is known that different hydrocarbons, for 
example aromatics and aliphatics, have different AH*. On the other hand, inert gases could 
be used as a test for the present model. 

Thus, according to our analysis, the converging values for entropy and enthalpy are the 
consequence of point solvation energy. This is an interesting conjecture. A rigorous statistical 
mechanical treatment of this problem seems possible, but has never been developed. In 
general, point solvation energy is dependent upon whether it is solid, or liquid, or gas, from 
which the solutes are transferred. A quantitative theory might be able to explain the small 
differences between the three SMPG plots. 10 (Note: after correction according to molarity 
scale, the intersection, AS*, for gas — > water dissolution is increased by about 60 J per 
degree per mole, see ref. 10). 

T* s and Tl 

It is natural to suspect that the unique convergence temperature, T*, for entropy of all 
three different groups of dissolution transferred from either solid, liquid, or gas, is due to 
some intrinsic properties of water. T* is the temperature at which AS(T) in Eq. (8) equals 
AS* given in Eq. (12). That is: 

(Jt + v^t = (14) 

i.e., 

By an approximated calculation, Ot/o"£ = 0.158/T. Hence: 

t>(t;)£(t;) = o.i58. 

This is consistent with the laboratory measurements T* = 383 K, £ = 0.363 and a = 
1.1 x 10~ 3 . 

Similarly for T£, from Eqs. (9) and (13) we have: 

a - To T - Ta^T = 0. (15) 

According to SPT, the first two terms cancel each other. The third term, as we have 
indicated, is indeed quite small. Thus unfortunately, Eq. (15) is buried in our various 
approximations and fails to provide an estimation for T^. A more accurate estimation is 
required to obtain the convergence temperature for AH. 
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Discussion 



With the physical insight provided by SPT, we now attempt to answer some key questions 
concerning the thermodynamics of hydrophobic solvation and the hydrophobic effect. 

1) What is the Hydrophobic Effect? 

Ever since Kauzmann's seminal paper on hydrophobic effect, 20 people have believed that 
hydrophobic effect is mostly due to the reorganization of hydrogen bonds among the solvent 
molecules around the solute, and the contribution of direct interaction between solvent and 
solute is rather minimal. So what is the relation between the hydrophobic effect and the 
solvation thermodynamics for hydrophobic solute in water? To address the question, one 
has to be precise about the meaning of "hydrophobic effect". There is an experimental 
(thermodynamic) side and there is a structural (theoretical) side of conventional wisdom on 
"hydrophobicity" . The experimental side is that non-polar solutes in water have very low 
solubility and the dissolution has large heat-capacity changes, in contrast to the dissolution 
in organic solvent. The structural side is that hydrogen bonding arrangement has been 
altered when an non-polar solute is dissolved in water. 

From our analysis, it seems that low solubility and large AC P in fact stem from two 
distinct sources. 21 While the large AC P is associated with rearrangement of water molecules, 
the low solubility is primarily due to the geometric properties of water molecules. This 
suggestion is consistent with our understanding of entropy-enthalpy compensation, which 
says that the ability of rearrangement of solvent should only have minor effect on solubility. 7 
The more dynamic aspect of water will be reflected only through quantities like entropy, 
enthalpy, and heat capacity. The reorganization process within solvent contributes to entropy 
change through heat capacity AC P , which in turn is related to the peculiar large temperature 
dependence of a, the thermal expansivity, of water. Hence it seems legitimate to identify 
AC P with the hydrophobic effect. 

However, could one simply identify the AC P term as hydrophobic free energy? This indeed 
is the central issue behind the work by Murphy et al. 10 It was, of course, well recognized 
that the AC P term is not uniquely determined until an appropriate reference temperature(s) 
is chosen. 22 Murphy et al. proposed the using of T* and as reference temperatures, 10 
and thereafter a AC P term was uniquely defined. However, as we have seen, the existence 
of convergence temperatures and their values are not the hallmark of the reorganization of 
associative solvent. We suspect that many other solvation processes might also have such 
convergent properties. The basis for the existence of convergence temperatures is Eq. (7), 
and it is clear that this equation is not unique for hydrophobic solvation (more discussion 
later). 

2) What is the Role of Hydrogen Bond? 
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An inevitable objection to our approach from many readers will be the complete ne- 
glect of hydrogen bond which has central importance in Kauzmann's structural model for 
hydrophobicity. 20 We would like to emphasize that we do accept the hydrogen bond in 
water as the structural base for hydrophobicity but we want to seek the specific ther- 
modynamic aspect or aspects of the hydrogen-bond structure which are responsible for 
the thermodynamics of hydrophobic solvation. The large temperature dependence of a no 
doubt is a manifestation of hydrogen bonding reorganization in water, the ultimate source of 
hydrophobicity. However, other aspects of the water molecules might also be relevant or even 
crucial; for example, the tetrahedral chemistry of hydrogen bonding, 20 or more generally 
the non-isotropic pair-wise interaction between two water molecules, 23 and the high ratio 
between physical volume and thermodynamic volume of water. 3 In our model, SPT indeed 
uses all these properties of water, though not explicitly. The fact is that water molecule in 
the SPT model has large thermodynamic volume, i.e., low packing density, but at the same 
time a small radius leads to high solvation number around a cavity. 3 This indicates that the 
solvent molecule is not isotropic, and there are preferences for these molecules to surround a 
cavity, which is exactly the Kauzamnn's argument! A crude analogy will be a wedge-shaped 
molecule, and that is quite consistent with water molecules. 

3) What is the basis of convergence temperature? As pointed out by Lee, 13 the presence of 
convergence temperature is due to some kind of linear free energy dependence on molecular 
substituents, i.e., Eq. (7). 

4) What determine the magnitudes of AS* and AH*? They are determined primarily by 
the thermodynamics of solvation of zero size point solutes. 

5) Why is there a large AC P ? The AC P of hydrophobic solvation stems from the peculiar 
large temperature dependence of thermal expansion coefficient, da/dT. 

6) How to Obtain Molecular Interaction Energy From the Thermodynamic Data? 

This is an age-old question. Twenty-five years ago, T.H. Benzinger proposed a new 
definition for enthalpy of chemical reaction. 24 His argument was that for a chemical reaction 
with non-zero AC p , there would be no unique heat of formation for the reaction. Since the 
true heat of formation is mechanical (athermal), Benzinger suggested to use AH(0) at zero 
Kelvin as the "true" heat of formation, and argued that the remains of free energy should 
be lumped into one term: 



AW(T) = AH(0) - AG(T). 
When taking into account the fact AS(0) = 0, he obtained: 
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In some sense, what Benzinger did was similar to what Murphy et al. did. They were both 
trying to divide the total free energy into a part with direct (mechanical) interaction and the 
rest part with surrounding effect due to thermodynamics. They both realized that AC P term 
was related to the latter since it characterizes the fluctuation in enthalpy due to thermal 
agitation. The crucial question is of course whether it is possible to find an appropriate 
reference temperature(s) based on purely thermodynamic analysis without any molecular 
model. Chan and Dill recently have extensively discussed this issue, and they concluded 
that without a molecular model, purely thermodynamic analysis would not provide much 
meaningful result. 12 

7. The Validity of Using SPT to Model Hydrophobic Solvation. 

Let's now reiterate the rationales for using SPT to model solvation in associative solvents. 
In SPT, the reorganization of solvent is considered implicitly through experimental data on £ 
as function of temperature. This approach is consistent with the assertion that the ultimate 
reason for reorganization in a solvent is its temperature dependence as a pure liquid. 7 The 
presence of hydrogen bonds between solvent molecules is manifested in the experimental 
data on water density and its temperature dependence. SPT, of course, completely neglects 
the soft interaction between solute and solvent. 

Our second defense is based on a recent analysis of how thermodynamic systems respond 
to small perturbations. 7 It has been shown that if we classify thermodynamic quantities 
by the orders of derivative of free energy, there is a relationship between thermodynamic of 
perturbation and the thermodynamics of unperturbed system. Since one can treat solvation 
as a perturbation, one only needs the thermodynamics of one order higher for pure solvent 
in order to calculate the thermodynamics of solvation. For example, ax of solvent gives 
AC P of solvation. All the structural changes in hydrogen bonding will be captured in these 
thermodynamic quantity of pure water, and our analysis made use of them. 

Finally, we would like to emphasize that we are not attempting to use SPT to model 
the properties of water; rather, we are merely using SPT to relate the thermodynamics of 
dissolution of an inert solute in water to that of pure water. In the past, many models for 
water which are based on multi-state of water conformation have been successful in providing 
calculations for thermal expansion coefficient, a, but have failed to deal with its temperature 
dependence. 25 This situation is completely in accord with our analysis. 
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Appendix: A possible AH* in gas to water dissolution? 

Let's proceed with the dissolution of a gas into water (g — > w) by first liquefying the gas 
(g — > 1) and then transferring the liquid into water (1 — > w): 

AHg^yj = AHg^l + AH . 

There seems to be a convergence temperature for dissolution enthalpy change from liquid 
— > water at about 20°C, and corresponding AH* is about few KJ/mol. 16 According to 
Trouton's rule, we can write 

AH g ^i = -88T b J/mol 

It is well known that the boiling temperature T b is scaled with molecular size, hence AH g ^ w 
and AHi^g should have similar converging AH* but at different temperatures. Mathemati- 
cally, if we have: 

AH l ^ w = AH* + b 1 (T)r 2 
where b x = when T = 20°C, and 

AH g ^ = b 2 (T)r 2 

then: 

AH g ^ w = AH* + [b^T) + b 2 (T)}r 2 
the convergence temperature is the value of T at which b\ + b 2 = 0. 
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